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(Model A). Jarzynski's equality is obtained through symmetries of the dynamical action in the 
path integral representation. We derive a set of exact identities that generalize the fluctuation- 
dissipation relations to non-stationary and far-from-equilibrium situations. These identities are 
prone to experimental verification. Furthermore, we show that a well-studied invariance of the 
Langevin equation under supersymmetry, which is known to be broken when the external potential 
\ is time-dependent, can be partially restored by adding to the action a term which is precisely 
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Ward-Takahashi identities. 
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I. INTRODUCTION 

During the last decade, many exact relations for non-equilibrium processes have been derived. The Jarzynski 
equality is one of these remarkable results: 
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(e-W) = e^ AF . (1) 



This relation implies that the statistical properties of the work performed on a system in contact with a heat reservoir 
at temperature kT = /3 _1 during a non- equilibrium process are related to the free energy difference AF between two 
equilibrium states of that system. This identity was derived originally using a Hamiltonian formulation [l| and was 
extended to systems obeying a Langevin equation or a discrete Markov equation [2|. This result was generalized by 
Crooks [1,01 j wn o showed that the identity ^ results from a remarkable relation between the probability Pp(W) of 
performing the quantity of work W in a given (forward) process and the probability Pr(—W) of performing — W in 
■ the reversed process, namely 

>: vf(w) _ 8(W _ AF) 

8: V R (-W)- e (2) 

00 ■ 

Jarzynski and Crooks' identities are now well established results (a review of the state of the art can be found for 
• ' example in 0]). These relations have been verified on exactly solvable models @ and by explicit calculations in kinetic 
theory of gases @, Q ■ These equalities have also been used in various single-molecule pulling experiments (9W1TI] to 
measure folding free energies (for a review of biophysical applications see e.g. [HI), and have been checked against 
analytical predictions on mesoscopic mechanical devices such as a torsion pendulum [l3j . Experimental verifications 
are delicate to carry out because the mathematical validity of Jarzynski's theorem is insured by rare events that occur 
with a probability that typically decreases exponentially with the system size [T3 |. 

Another type of physical systems where large fluctuations are expected to occur are extended statistical models 
in the vicinity of a phase transition. Such systems are often modeled by a continuous space-time description with a 
local coarse-grained order parameter, 4>(x), which minimizes a Ginzburg-Landau type free energy. The equilibrium 
properties of these models have been thoroughly studied, in particular using renormalization group techniques [T^-[l3 ■ 
Besides, in the vicinity of a critical point, the dynamic properties also display anomalous behaviour. If one assumes 
that universality remains valid 16], it is natural to construct and investigate the simplest dynamical models, with 
a given static behaviour, which respect some physical constraints such as symmetries and conservation laws. The 
coarse-grained order parameter of a microscopic system is then represented by a space and time dependent field <j>(x, t) 
that evolves according to an effective stochastic differential equation. Different possible types of evolution equations 
have been classified (see e.g. the review paper of Halperin and Hohenberg [18|). 
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In the present work, we derive nonequilibrium work relations for a field 4>{x, t) that follows the simplest dynamics 
in Halperin and Hohenberg's scheme: Model A dynamics describes the kinetic Ising model with non-conserved order 
parameter. We represent the stochastic evolution as a path integral weighted by a dynamical action [13] ■ Our method 
is closely related to the one used in [2(], [2l[ to study the case of a Langevin equation for a O-dimensional scalar 
coordinate that depends only on time. Adding a spatial dependence allows us to use the powerful response-field 
formalism that was developed in [ljj. The work relations are then derived from elementary invariance properties 
of the path integral measure under changes of variables that affect simultaneously the original field <p(x, t) and its 
conjugate response- field 4>{x, t). From the work relations, we derive correlator identities that generalize the equilibrium 
fluctuation-dissipation relations for situations that can be arbitrarily far from equilibrium. 

One advantage of introducing space and time varying fields is to extend the possible symmetries of the system and to 
consider transformations that can mix space and time. In the present context, this can be achieved by introducing two 
conjugate auxiliary Grassmann fields, c(x, t) and c(x, t). The new dynamical action, which now depends on four fields 
4>, (f>, c and c, exhibits a larger invariance which is a manifestation of a hidden supersymmetric invariance [22I [23| . This 
property is true at thermodynamic equilibrium and it is known that the equilibrium fluctuation-dissipation theorem as 
well as the Onsager reciprocity relations can be derived from it (24j-[26| (see also [l7j )■ In other words, supersymmetry is 
a fundamental invariance property of the full dynamical action that embodies the principle of microscopic reversibility. 
For a system out of equilibrium (for example a system subject to a time-dependent external drive), supersymmetric 
invariance is broken. This leads to a violation of the fluctuation-dissipation theorem or, equivalently, to the occurrence 
of corrective terms in the formulation of this theorem: this fact was clearly recognized in [271 281 ]. Here, we remark 
that weighing the expectation values by the Jarzynski term e~^ w amounts to modifying the dynamical action of the 
model and we show that the modified action exhibits an invariance under a specified supersymmetric transformation. 
This invariance manifests itself as correlators identities known as the the Ward-Takahashi identities (a field-theoretic 
counterpart to Noether's theorem). Finally, we prove that the nonequilibrium work relation can be deduced from 
the Ward-Takahashi identity that encodes the underlying supersymmetry. Therefore, supersymmetric invariance of 
stochastic evolution equation is a fundamental property that embodies equilibrium relations (Onsager reciprocity, 
fluctuation-dissipation theorem) as well as nonequilibrium work identities. 

The outline of this work is as follows. In Section HH we define the model, use the response-field formalism to 
derive the work relations for a space-time dependent field, and obtain a fluctuation-dissipation relation, valid far 
from equilibrium, which can be shown to be mathematically equivalent to the Jarzynski relation. In Section IIII1 
we use the formalism of Grassmann fields for the equilibrium case and define precisely the various supersymmetric 
transformations that leave the dynamical action invariant. We then show that supersymmetry, which is broken when 
the system is out of equilibrium, is restored by modifying suitably the action and we prove that Jarzynski's equation 
can be viewed as a consequence of the Ward-Takahashi identity that encodes the restored invariance. Concluding 
remarks are given in Section HVl Technical details are given in the appendices. In particular, the superfield formalism 
and the derivation of the Ward-Takahashi identities are recalled in Appendices |B] and [Cj 



In this section, we derive the field-theoretic version of the non-equilibrium work relations for a system that obeys 
a time-dependent Ginzburg-Landau equation. The dynamics considered will be purely relaxational and we shall 
focus on the most elementary case with no conservation laws, described by Model A dynamics. We shall express the 
Probability Distribution Function of the field at a given time as a path-integral. The Jarzynski and Crooks relations 
will be obtained via this path integral formalism. 



II. STOCHASTIC EVOLUTION OF A SCALAR FIELD 



A. Model A dynamics 



We consider a scalar field <fr(x,t) that evolves in a d-dimensional space according to Model A dynamics [HI, given 
by the following stochastic equation of motion: 




r 



6U[<j>(x,t),t] 
5<fi(x, t) 



(3) 



where the dynamics is governed by the time-dependent potential 




(4) 
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h(x,t) being an external applied field. The time-independent part of the potential assumes the familiar Ginzburg- 
Landau form and is given by 

F GL [<j>] = J A d x{ l -r a cf + ||W| 2 + y0 4 } ■ (5) 

The fluctuating driving field £(x, i) is assumed to be a Gaussian white noise of zero mean value and of correlations 
given by 

(C(x, t)C(x',t')) = 2^6(t - t')8 d (x - x') , (6) 

where f3 — (fcT) -1 is the inverse temperature. The fact that the auto-correlation of the noise satisfies Einstein's 
fluctuation-dissipation relation ensures that the dynamics is microscopically reversible and obeys detailed-balance 
[30( | . The Langevin equation ^ for model A thus reads 

^(x, t) = -T f{<j>{x, t),t)+ C(x, t) , (7) 
where, for later convenience, we have defined 

fMx, t), t) = SU ^^ t] = -V 2 0(x, t) + r <f>(x, t) + uo^ix, t) - h(x, t) . (8) 

The dynamics of the order parameter can also be described in terms of the Probability Distribution Function (PDF) 
V(4>i \4>o) of observing the field <f>{x,t) = 4>\{x) at time t = tf, knowing that the initial field is 4>o{x) at time t = 0. 
By definition, this PDF is given by 

V{4> 1 \<i> ) = (5{ct>{x,t f )-<t>i{x))) ) (9) 

where the expectation value is taken over all possible realizations of the noise ((x, t) between the initial and the final 
times (the initial value of the field 4>o(x) being fixed). Substituting the Gaussian measure for the noise, this expression 
becomes 

V(.<h\<h) = I V((x,t) e~^ fddxdtc2 S(<b{x,t f )-Mx)) with <p(x,0) = cj> (x) . (fO) 



This expression is nothing but a formal path-integral solution of the functional Fokker-Planck equation associated 
with the Langevin dynamics (J3J): 



f-„/^(^>p + if). 



When the external field is constant in time h(x,t) = h(x), this Fokker-Planck equation has a stationary solution, 
which is the equilibrium Gibbs-Boltzmann distribution: 

P -0U[<i>] r 
P *M = -T^nrTT with Z iP> h)= V<j> e-ww . (12) 



eq ^> z\p,h] 

Finally, we recall that the equilibrium free-energy F[/3, h] is defined by 



F[l3,h} = --log Z[0,h}. (13) 



B. Dynamic action for the Probability Distribution 



The probability V(4>i\4>o) of observing the field 4>\{x) at time tf starting from <f>o(x) at time t = is given by 
equation (|10p . We now rewrite the path integral in terms of the variable <p(x,t) using the response-field formalism 
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of Martin-Siggia-Rose, de Dominicis-Peliti and Janssen [l9|. We start with the following identity (that can be found 
e.g. in [13) 

r r<p(.x,t f )=cl> 1 (x) / 5U \ 

1= VMx) Vct>(x,t) S[<t>(x,t)+r —-{(x,t) |detM| 

J J<f>(x,0)=<j>o{x) \ °<P J 

V4>{x,t) 5 U(x,t)+r ^-COM)) |detM|. (14) 

:x,0)=4>o(x) \ 0<P J 

Note that in the second equality, the configuration (j){x,tf) of the field at time tf, appears as an integration variable 
(i.e. < t < tf). The linear operator M is defined as 

SC(x,t) d df{4>{x,t),t) 

M= 6^J) = di + Tn W • (15) 

The determinant of this operator can be written as 

detM = exp{Tr(logM)} = e 2 J ^ , (16) 

where the last equation is found by discretizing the operator M and using the Stratonovich convention [I?], Hd| ■ We 
substitute this expression in the identity (TUT) and introduce the response field (j>{x,t) that allows us to rewrite the 
functional Dirac distribution 8{) as an exponential. Thus, we obtain 

1= f V<f>(x,t)T>0(x,t)\detM\ e -/^W+r ^-£} where < t < tf ( 17 ) 

We now insert this identity (fT7|) in equation f| 10[) . perform the Gaussian integral over the noise variable C( x it) an d 
substitute the expression (|16l) for the Jacobian of M. Finally, the following expression for the PDF is obtained: 

r<t>{x,t s )=<pi{x) 
J 4>{xfi)=4>a{x) 

The PDF is thus expressed as a path-integral over the order parameter 4>(x,t), with an effective dynamical action £ 
given by 



SU_ TqS 2 U 
T ' ^"^ )_ T 



E(0,0^) = r o 0(^ + — . (19) 



The non-equilibrium identities will arise from invariance properties of the path- integral with action E. 



C. Non-Equilibrium Correlations Identities 



We consider the case where the applied field varies with time according to a well-defined protocol: for t < 0, we 
have h{x 1 0) = ho(x) and the system is in its stationary state; for t > 0, the external field varies with time and reaches 
its final value h f (x) after a finite time t f and remains constant for t >tf. The values of the potential IA for t < and 
t > tf are denoted by Uq and U\, respectively. 

Let 0[4>] be a functional that depends on the values of the field 4>{x,t) for < t < tf. The average of 0[cj)] with 
respect to the stationary initial ensemble and the stochastic evolution between times and tf is given by the path 
integral 

(0) = \- ( V ( f) {x)V^ 1 {x)e' l3Ua ^ ' ^ M ) p^e^ ddxdtsw <^ ) C>[0] 



Z 



x,0)=<Mx) 



„-/3W o [0o] 

= / V, .,.,-)/ Vfo(x) ^ <0i|g>|0o) (20) 



where we have defined 



(0i|O|0o)= [**'*' V4>(x,t)V^(x,t)e-^ dxdt ^^ 0[<f>] . (21) 
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Under a change of the integration variable in equation (|2"01 , 



the path integral measure is invariant but not the action S which varies as 



(23) 

o<p[x, t) 



Noticing that 



we obtain 



■ 8U[(j>(x,t),t] _ dU[(j){x,t),t} dU[<j>(x,t),t) 
5(f>(x,t) ~ dt dt 



(24) 



ddx r*t dt ^suypM = 

Jo 0(j)(x,t) 



The last term in this equation represents Jarzynski's work, defined by 



rtf oil 

WjM = / dt — = - I d d xdt h{x,t)ct>(x,t) , (26) 



the last equality being a consequence of equation Q . The change of sign of the time derivative <f> in equation 
now compensated by the change of variables in the path integral 

(<f>{x, t),$(x, t)) -> (<Kx, t f - t), 4>{x, t f - t)) . (27) 

This time-reversal transformation leaves the functional measure invariant and restores E to its original form but with 
a time-reversed protocol for the external applied field h(x, t) — > h(x, tf —t). Performing the above change of variables 
(|2"21 and (j2"7| in equation ([2~Tj) and using equations (|2"5|) and (j2"5)) . we find, recalling that the work Wj is odd under 
time-reversal, 

{^O^o) = e^ Uo ^- Ul ^^(t>i\e-^ Wj 6\(j) ) R . (28) 

On the right hand side, the subscript R on the expectation value denotes a time-reversed protocol. The notation with 
a hat "over an operator denotes the time-reversed operator, more precisely: 

6{4>} = 0{<t>{x,t f -t)}. (29) 

Inserting this identity in equations ((201 - 121]) allows us to derive the following general relation: 

(O) = / VM^VM^e-^HMe-^Ol^n = ?±(6e-^) R = e^ AF (6e-^) R , (30) 

where AF is the free energy difference between the final and the initial states. Finally, redefining O as Oe~@ Wj , we 
deduce that 

(Oe-Wj) =e-? AF (d) R . (31) 

When O — 1, we obtain Jarzynski's theorem 

( e -W) = e -/JAF_ (32) 

Taking O = e {l3 -^ Wj , where A is an arbitrary real parameter, we derive the following symmetry property 

(e- xw ')=e- 0AF (e^-V w ') R . (33) 
The Laplace transform of this equation leads to Crooks relation ([2j in its usual form 0, 0] : 

Vp{W) _ o/3 ( W _AF) 



V R (-W) 



e PKw-*n (34) 
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where Vf and Vr represent the probability distribution functions of the work for the forward and the reverse processes, 
respectively. We emphasize that the proof of Crooks and Jarzynski identities is based on invariance properties of the 
path integral and does not involve any a priori thermodynamic definition of heat and work. The expression (|26|) for 
the Jarzynski work appears here as a natural outcome of this invariance. It is important to notice that time-reversal 
is crucial to obtain the general identity (|31l) and Crooks' theorem. However, it is known that Jarzynski's identity can 
be proved without assuming time- reversal invariance 0] . 

We emphasize that Jarzynski's identity is valid only under carefully defined boundary conditions: (i) the system is 
at thermal equilibrium at time t = 0; (ii) During the finite time interval < t < tf, the system is subject to an 
external protocol that drives it away from equilibrium; (iii) After the finite time tf, all time-dependent parameters 
are frozen: these fixed parameters define a new state of thermodynamic equilibrium towards which the system relaxes 
after an infinite amount of time. According to this scheme, all path integrals must range over the finite interval of 
time < t < tf and the expectation value of the operator O, defined in (f^Uj) . has to be taken with respect to the 
Boltzmann distribution at t — and the uniform distribution at the final time tf. However, these stringent boundary 
conditions necessary for Jarzynski's identity to be valid, allow us to embed naturally all the path integrals over the 
infinite range of time — oo < t < +oo by using the following properties of the probability distribution: 

1 -W0o] = lim p(0 o |^ t ) ; (35) 

ZjQ T— > — OO 

1 = /" Vc/>(x,T)P((f) T \(f>t f ) for any r > t f . (36) 

The first property assumes ergodicity (i.e. the Gibbs-Boltzmann distribution is reached at time t = by starting from 
any initial condition at t = — oo). The second equality simply results from normalization. In terms of path integrals, 
the first expression becomes 



J_ e -/3W o [0o] 



z 



0(x,O)=0 o (x) 

V<p(x, r)£>0(x, r) e" / d xdts W.^) for - oo < r < , (37) 

'(x, — OO) — (/j — oofx) 



where the condition at t = — oo is taken to be an arbitrary value 4>-oo (or more generally a distribution of values, 
normalized to 1). Similarly, equation (1361) is rewritten as 

1=1 T>ct>(x,T)V4>(x,T)e- / dd * dts (^) where t f < r < oo . (38) 

J 4>{x,t 1 )=4> 1 {x) 

We now consider an operator 0[<f>] that differs from a constant only for < t < tf (i.e. the operator 0[(j>] depends on 
the values taken by </> only over the finite range of time < t < tf). Using the relations ([37| and (l38l) . the expectation 
value of 0[(j>], defined in equation (flHI)) . can be expressed as 

(O) = I V(t>( X ,T)Vj)(x,T) e -/d d *dtE(^,S> foj . — OO < T < OO , (39) 



and where the space-time fields 4>(x, t), 4>(x, t) are integrated over an infinite range of time and over the whole space. 
The only restriction on this path integral is that the initial condition at t = — oo is fixed (or more generally, the values 
at t = — oo are sampled from a normalized distribution). We note that the Jarzynski term e~^ Wj is equal to 1 outside 
the interval <t <tf and therefore we can also write 

/ p -PWj\ _ f TtA-DA P -Sd d xdt-£(<f,,j,,$) p -pWj 



We emphasize that the boundary conditions at finite time, as well as the average over the Boltzmann factor, have 
been eliminated and all path integrals are now evaluated over the full time line — oo < r < +00. 



D. A Non-Equilibrium Fluctuation-Dissipation Relation 

The identity (|31j) . which is at the core of the work fluctuation relations, is valid for any choice of the external field 
protocol. The free energy variation is a function only of the extremal values of the applied field at to = and t = tf 
and is independent of the values at intermediate times. Therefore, performing functional derivatives of the Jarzynski 
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identity (1321) with respect to h(x,t) at an intermediate time to < t < tf, and at position x, results in new identities. 
For example, we have 

More generally, the n-th functional derivative of equation (|32[) at intermediate times t\, . . .t n and positions x\, . . . x n , 
gives the identity 



'e-W'Hifouti) ~ fkxuti))) = 0. (42) 
i o 



Similarly, the functional derivative of equation (I31[) leads to 

($(x,t) - £^(x,t))Oe-^) = e-^ AF ($(x ) t)d) R . (43) 
t o 

In particular, equation (|4"Tj) follows by choosing = = 1 and taking into account that ((f)) = 0. [Indeed, if we take 
the functional derivative of equation (|2"0)l with respect to h(x,t) for O = 1, we obtain (0) = (jP + j^) = 0.] 
For the special case 0[cj)] — cj)(x',t'), equation (|4"3"]l leads to: 



(4>{x, W(x>,t')) - ^-{j>(x, t)<f>(x', t')e-^) = e-? AF {$( Xl t)4>(x', t')) R , (44) 
J- o 

where 4> is obtained from 4> by time reversal as defined in equation (|29[) . The terms proportional to the response field <f> 
in the correlators can be generated as follows: we consider a small perturbation h\(x, t) that drives the system out of 
the fixed protocol h(x, t). However, we keep the definition of the Jarzynski work unchanged so that the perturbing field 
h\(x,t) is not included in Wj. The field hi couples to <j) in the action S: therefore performing functional derivatives 
with respect to h\ amounts to inserting the field <fi inside correlation functions. The previous equation can thus be 
rewritten as 

r dhi(x,t) h!=o 5hi(x,t) 'u=o 

We emphasize that Wj has to be measured for the fixed protocol h(x,t). In this form, equation (|4"5l) appears as 
an exact generalization of the fluctuation-dissipation relation (FDR). The equilibrium FDR [3ll - l33j is retrieved by 
setting Wj and AF to 0. This indeed corresponds to a system prepared in an equilibrium state, which is not subject 
to any macroscopic protocol (i.e. h(x 1 1) — 0) but which is driven slightly out of equilibrium by the small perturbation 
h!(x,t): 

f(<t>(x>,t>)<j>(x,t)) = (ct>(x',mx,t)) - (4>(x',t')<f>(x,t)) = 5 f h {X '^ - S^'/Vx . (46) 
I o dhi(x,t) ' 



^i=o 5hi(x,t) 



hi=0 



The fact that the equilibrium fluctuation-dissipation relation can be deduced by diffentiation from Jarzynski's identity 
(or equivalently from Crooks' relations) has been understood by various authors (see in particular the works of 
R. Chetrite et al. [34T - [37| ). This technique can also be used to find analogs of the FDR at higher orders [Hj|. 
Generalizations to non-equilibrium stationary states (NESS) have been also proposed 13a. I39l - l41| . e.g., starting from 
the Hatano-Sasa relations which are the counterpart of Jarzynski's identity for a NESS [42|]. We emphasize that the 
identity obtained in equation ((45]) belongs to a different class. We do not consider a linear perturbation near a state of 
thermodynamic equilibrium, or near a NESS. Rather, we first apply, as in Jarzynski's scheme, a protocol to a system 
initially in thermodynamic equilibrium (that can be driven as far from equilibrium as wished) and then, we apply 
linear perturbations around this fixed protocol: this leads to a new fluctuation-dissipation theorem that relates out 
of equilibrium and nonstationary response functions to nonequilibrium and nonstationary correlation functions. The 
insertion of the Jarzynski factor e~@ Wj inside the correlators leads to formulae which are valid far from equilibrium and 
look very similar to equilibrium relations. The relation (|45[) could be verified in single molecule pulling experiments 
where the protocol corresponds to the pulling force F(t) and <fi does not depend on space. Then, all the quantities 
that appear in this relation are susceptible to experimental measurements by adding a small perturbation 5F(t) to 
the fixed protocol F(t). 

The correlator identity (|4ip was obtained as a consequence of Jarzynski's equality (f3"2"j) by taking its first derivative. 
Conversely, we show in Appendix El that equation (|41[) implies Jarzynski's equality (|32[) and is therefore equivalent 
to it. This converse property will be used in the next section to show that the work relation can be extracted from a 
hidden supersymmetric invariance of the dynamical action. 
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III. SUPERSYMMETRY AND NONEQUILIBRIUM WORK RELATIONS 

Identities between correlators such as equations pT1) - ([4"5)) suggest the existence of an underlying continuous sym- 
metry of the system. Indeed, it was recognized in the late seventies that the Langevin equation possesses a hidden 
invariance under supersymmetric transformations. This property was first discussed in the context of dimensional 
reduction [22| and then used to derive convenient forms for diagrammatic expansion techniques (23j that were used 
to study critical dynamics of relaxational models (43l - l45l |. Supersymmetry was also used to st udy mesoscopic quan- 
tum systems (see e.g. in (4(| H?})- This supersymmetry became an efficient tool to study [48Tl5lj | the properties 
of Fokker-Planck and associated Schrodinger operators (see e.g. [H, HU). It was also realized that the equilibrium 
fluctuation-dissipation relation and the Onsager reciprocity relations could be derived from this invariance |24| - [26l l44j . 
Conversely, nonequilibrium situations were found to correspond to supersymmetry breaking and corrections to the 
classic equilibrium relations could formally be calculated (27], [28| . In this section, we extend this investigation further 
by showing that although supersymmetry is broken under nonequilibrium situations, it is partially recovered by adding 
to the dynamical action a term, which precisely corresponds to Jarzynski's work. This restored invariance leads to 
Ward-Takahashi identities amongst correlation functions. Jarzynski's relation results from these identities. Thus, it 
is emphasized here that supersymmetry is not imposed but it is rather a result, namely an underlining symmetry of 
this system. 



A. Supersymmetric Invariance for the time-independent Model A 

First, we consider the case where the external field does not depend on time: the Langevin equation Q has then 
a well-known supersymmetric invariance. We shall review the formalism that allows to make this invariance explicit, 
write the Ward-Takahashi and derive the equilibrium fluctuation-dissipation relation following 24 1. 



1. The Supersymmetric Action and its Invariance Properties 

To uncover this hidden symmetry, we introduce in addition to the original field 4>{x, t) and the response field 4 > {x 1 1), 
two auxiliary anti-commuting Grassmann fields c(x. t) and c(x, t) that allow us to express the Jacobian of M, defined 
in equation flT5l) . as a functional integral [13, [2J, |25[ . These fields c{x,t) and c(x,t) can be viewed as hidden classical 
fcrmionic fields that ensure the volume conservation constraints: they allow to enforce this conservation property at 
a dynamical level. Inserting the following identity [17| 

detM= / VcVce cms = / VcDce { ^ +T °^ Y \ (47) 



we observe that the PDF can be rewritten as 

n(j>(x,ti)=(l>i(x) _ 

P(0i|0o)= / V<j)V(f)VcVc e~ f d xdts (<^.c,c) , (48) 

J 4>(xfl)=4> (p) 

where the effective Lagrangian S, which is now a function of the Grassmann variables as well, is given by 

E(**^-r c) «i + *-f)-4 + r„gk. (49, 

The action S exhibits two important invariances under infinitesimal transformations that mix ordinary fields with 
Grassmann fields. We shall now describe them by specifying how each field varies under these transformations. In the 
Appendix [B] we shall use a more elegant presentation in which the four fields (</>, <j), c, c) appear to be the components 
of a unique superfield <&; also, in this language, the dynamical action S will take a more compact form and the 
infinitesimal transformations that leave it invariant will have a simple interpretation. 



• Invariance under BRST1 Transformation: Consider e to be a time-independent infinitesimal Grassmann 
field. We consider the following transformation (that we call BRST1): 



5<fi(x,t) = —c(x,t)e, Sc(x,t) = , 

5c(x,t) = 4>(x,t)e, S(f>(x,t) = 0. 



(50) 
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We note that the square of this transformation vanishes. If we calculate the variation of X under the transforma- 
tion (jSTJ)) . we obtain using equation (|39"]k 



= 0(4 + r vtj 8<f>) -Sc(— + T — )c - T 0J -j5<l> cc = -0(5e + ^ —ce) - 0e(t+T Q —c) = . (51) 



This expression vanishes identically because of algebraic anti-commutation rules. We note that we do not need to 
suppose that the potential U is time-independent. 

• Invariance under BRST2 transformation: The transformation BRST2 mixes the different fields as follows : 

5<f>(x,t) — c(x,t)e, 5c(x,t)=0, 

5<j>(x,t) = ^-c(z,t)e, 6c(x,t) = \$(x,t) - ^<j>(x,t)je, (52) 

e being a time-independent infinitesimal Grassmann field. Here again, the square of the transformation (|52j) vanishes. 
If we calculate the variation of £ under the transformation (1521) , wc obtain 



„ d f//3 . -\ 1 JU. J 2 U . d (fj3 . JU A 1 Q 5 2 U 



If the potential U is independent of time the last term vanishes and S is invariant under the transformation (|52p only 
up to a total time-derivative term that produces boundary contributions to the total action. This does not affect the 
dynamics if boundary terms vanish or if the time-integral is defined from — oo to +oo. 

The invariance of the dynamical action £ under both transformations BRST1 (|5TJ1) and BRST2 (|5"2"|) is what makes 
the time- independent dissipative Langevin equation supersymmetric [17j . This supersymmetric property reflects the 
time reversal invariance of Model A in the absence of an external field and allows to prove the fluctuation-dissipation 
theorem [24], [2f| as will be recalled below. 

It is important to emphasize here that in the general case the use of time reversal invariance in a path integral of the 
type in Eq. (|48|) . as mentioned above, extra precaution is needed. Whereas the topology seen in Ea. (|48l) is trivial, the 
general case of non-trivial topology is subtle and time reversal invariance arguments must take this into account. The 
case of non-trivial topologies was studied in details in ref. [54] where large deviations and global topological currents 
were discussed in the case of an Heisenberg spin-chain with a Wess-Zumino type term. Other work on non-trivial 
topologies can be found in ref. 55] 

2. Ward-Takahashi identities and the equilibrium fluctuation- dissipation relation 
Introducing a four-component source (H, H , L, L) , we define the generating function 



Z(H, H, L,L) = J V(f>V<j)VcVcexp yj c\ d xc\t (-£(</>, 4>, c, c) + Hcf> + H<j) + Lc + Lc) 

In the Appendix, we rederive the following two Ward identities that result from the invariance of the action under 
the transformations ([5"Pf and (j5"2"j) . The first Ward-Takahashi identity corresponding to invariance under (|5"P|) is 



d d xdt [h S 4-L^-)=0. (55) 



5L SH 

The second Ward-Takahashi identity corresponding to the transformation ([52]) is given by 



f lrf , ( (3 d SZ T fSZ (3d8Z\ -SZ\ 

J dxdt {i H a-tJL +L {m-idtm) +H TL) = 



h Jl)=°- (56) 

We now apply 5 2 /SL(x\ t')5H(x, t) to the first Ward-Takahashi identity (1531) and then put all the sources H, H, L, 
and L to 0. This leads to 

S ' Z S ' Z =0. (57, 



5L(x', t')5L(x, t) SH(x', t')SH{x, t) 
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This equation implies, using (|C8|) and (|C9|) . the following identity between correlation functions: 

(c(x',t')c(x,t)) = ^(x',t')^(x,t)) . (58) 

Similarly, applying 5 2 /SH(x' ,t')6L(x,t) to the second Ward-Takahashi identity (l56l) . and putting all the sources to 
zero, we obtain 



8H(x',t')5H(x,i) T SH(x',t') {dt SH(x,t) J 8L{x,t)SL(x',f) 
This identity implies that 

<0(x',O#M)> " f&ix'xAx^)) + (c(x,t)c(x',l/)) =0. (60) 

r dt 



In order to eliminate the correlations between the Grassmann variables, we combine equation (|57p with equation (|59[) 
(or equivalently equation (|58p with equation (1601) ) and use the fact that L and L anti-commute (or equivalently that 
c and c anti-commute). This leads us to: 

f o (cf>(x',t')^x,t)) = (cf>(x',t>)t(x,t)) - <<H*',i')#M)> = (<M*',0> - aff(x',*') <#M)) ' (61) 

Recalling that 5 (4>) /SH is a response function, we observe that this equation is nothing but the Fluctuation-Dissipation 
relation (see (1461) ). Usually, the FDR is derived by invoking invariance under time-reversal which implies detailed 
balance. Here, it is the invariance under supersymmetry that plays the role of time-reversal invariance. 

B. Model A with a time dependent potential 

We now study the case where the potential U[4>(x,t),t] that appears in Model A ((3]) depends explicitly on time, 
and show how properties related to supersymmetry can still be used in this nonequilibrium situation. 

1. Breakdown of the invariance for a time- dependent potential and its restoration by adding the Jarzynski term 

When the potential I4[<p(x, t), t] depends on time (by adding for example a time-dependent external field), the action 
X is no more invariant under supersymmetry. More precisely, we observed that invariance under (|50l) does remain 
valid even when U is a function of time and therefore the first Ward identity ([53]) is still satisfied. However, invariance 
under (|52[) is broken and according to equation (I53[) . we find the variation of S(0, 4>, c, c) to be 

= d i-4 t {%)^), (62) 

with the total derivative term 

Therefore <5E is not a total derivative and invariance under (I52[) is not true anymore for a time-dependent potential. 
In particular, the second Ward identity (|56[) . which was crucial for the proof of the Fluctuation Dissipation relation, 
is no more satisfied. 

However, we note that the last term in equation (1621) . which breaks the invariance under (|52p . can be rewritten as: 

and can be interpreted as the variation of a function. Therefore, the modified action S j , defined as 

£ J = £ + /^ (65) 
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and obtained by adding the Jarzynski work (|26p to the initial action, is now invariant under the supersymmetric 
BRST2-transformation (|52j) because its variation is given by a total derivative term: 

= d i m 

However, we emphasize that Sj is no more invariant under BRST1 (1501) although S was invariant. We have thus 
restored BRST2-invariance at the expense ol BRST1. Therefore, in the time dependent case, neither the original 
action £ nor the modified action Sj are supersymmetric (l7j . They only exhibit partial invariances by either BRST1 
(in the case of X) or BRST2 (in the case of Sj). We shall now see that BRST2-invariance is required to derive 
nonequilibrium work identities. 

For a time-dependent external field h(x,t), the compensating term in equation (|64[) is given by —fih(x,t)(f). The 
boundary terms at t = ±00 are, conventionally, assumed to vanish. Therefore, the invariance of the dynamical 
action under (|52p breaks down when the potential is time-dependent but is restored by adding Jarzynski term. This 
observation allows us to use the Ward-Takahashi identity (|56p that results from this invariance. We shall prove that 
this Ward-Takahashi identity leads to the nonequilibrium work relations. 



2. Work relations from super symmetry 



We now show that the invariance of Sj implies the correlator identities (|42j) . We first remark that in the above 
proofs of supersymmetric invariance it was noted that the boundary terms (a total time-derivative contribution) are 
harmless if the integration range of the path integral is from —00 to +00. 

We consider an operator 0[4>] that differs from a constant only for < t < tf. Then, as shown in equations ([55]) 
and (|40p . the expectation value (Oe~^ Wj ) can be rewritten over an infinite time range. Using Grassmann variables, 
we have 

(Oe-' 3Wj ) = J V(f>V4>VcDc e -f^ d ^ti:(<P^c,c) e -Dw. 10 ^ _ J v^VcDc e~ S d^dtEjC^.c.c)^ ; ^ 

where all the space-time fields (f>(x,t), <j)(x,t), c(x,t) and c{x,t) are integrated over the time interval — 00 to 00. 
In the last equation, we have combined the action S with the Jarzynski work Wj to get the modified action Sj, 
defined in (|65p . As shown in the modified action Sj is invariant under the transformation ([52")) . Therefore, the 
following generating function Zj(H, H, L, L), built from the modified action Sj, satisfies the second Ward-Takahashi 
identity f5B]): 

Zj{H, H, L,L) = J V<j)V4>VcVcexp (J d d xdt 4>, c, c) + H(j> + H4> + Lc + Lc)^j . (68) 

We apply the following operator to the Ward-Takahashi identity ([56]) satisfied by Zj, 

6 TT ( 5 13 d 6 \ 



SL(x,t) jJ^\5H(xi,ti) r dti 5H(xi,ti) 
and set the source fields H, H, L, L to zero. For n = 1, we find 

^(x,t)-fj(x,t))e-? w ^=0 (70) 

More generally, for n > 1, we have 

Ufa - 7^0i)(02 - ffo) ...{4>n- ^-k>- pWj ) = , (71) 
\ J- to t / 

where <j)\ = 4>(x\,t{) etc... These two relations are identical to equations (|4"Tj) and (|4^1) . respectively. In Appendix \K[ 
we show that these relations are equivalent to Jarzynski's identity. This concludes the proof that Jarzynski's relation 
can be obtained as a consequence of a Ward-Takahashi identity that itself results from supersymmetric invariance. 
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IV. CONCLUSION 



We have used field-theoretic methods to derive nonequilibrium work identities for a space-time field driven by a 
non-linear stochastic equation (Model A). We have obtained a generalization of the fluctuation-dissipation relation 
that remains valid far from equilibrium and that characterizes the response of a system to infinitesimal perturbations 
around a given protocol. The introduction of auxiliary fermionic fields has allowed us to explore general symmetries of 
the dynamical action. In particular, it is well-known that the time independent Langevin equation exhibits a hidden 
supersymmetric invariance [22I [23[ that is known to imply the classic fluctuation-dissipation theorem and Onsager's 
relations [I?], HH [H| ■ However, this invariance breaks down when the potential varies according to a time-dependent 
protocol and drives the system out of equilibrium. In this work, we have shown that the invariance of the effective 
action under supersymmetric transformation is restored by adding to the action a counter-term which is precisely the 
Jarzynski work Wj. Furthermore, we proved that the associated supersymmetric Ward-Takahashi identity implies 
Jarzynski's theorem. Hence, supersymmetry enforces the exactness of the adiabatic limit even for processes that 
have a finite duration and that can bring the system arbitrarily far from equilibrium. In other words, weighing all 
averages with the Jarzynski work (which amounts to modifying the dynamical action by adding to it the Jarzynski 
work, as in (I65p ) restores one of the fundamental symmetries valid in equilibrium. Thanks to this invariance, many 
properties of the system are effectively the same as if it were at equilibrium (although it is neither in equilibrium 
nor in a stationary state). The idea of considering weighed averages, or equivalently modified path-measure (as 
was emphasized by Jarzynski himself in his early works [2j), allows to preserve certain crucial symmetries and has 
striking consequences in the present context. We believe that a similar arguments should apply in many different 
fields: in particular, supersymmetry exists in classical Hamiltonian systems [2r| for which Jarzynski's equality was 
initially derived, and can be applied to prove the the fluctuation theorem for stochastic dynamics [561 ] . Besides, the 
response-field method that we have used here can be extended to multi-component fields, to other stochastic models 
with conserved order parameter and also to systems with colored noise (57j . 
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Appendix A: Proof of the equivalence between equation (I41[) and Jarzynski's relation 

In this Appendix, we prove that the correlator identity (|4ip implies Jarzynski's relation (f3"2"j) and is therefore 
equivalent to it. First, we modify the applied external field h(x,t) by considering h(x,at) for any a > 0. We then 
evaluate the average value (e~" Wj )(a) using expression PU|) . From (|19l) and (|26p . we observe that the external field 
h(x, at) appears only in the following two terms: To4>(x, t)h(x, at) + /3ah(x, at)<f>(x, t) (note that the Jacobian term, 
5 2 U /Scf) 2 does not contain h). Thus, we have 



d{e 



-0W.J 



- /o^^W--^ /d^^t^^ + ^-j + .^-M,,.)) . (Al) 



da 

Integrating by parts the last term with respect to time leads to 
da 



= I p^ e -/3W J -/d d xdts(0,0,0) J d d xdt(T Q t$(x,t)h(x,at) - /3th(x,at)<j>(x,t)^ (A2) 
T>cl>V0e-P Wj -f ddxdts W'i'® f d d xdtth(x,at) (v ${x,t) - 04>(x,tj\ . (A3) 



Note that the boundary terms that result from integration by parts vanish because we are integrating for t G (— 00, +00) 
and because h = outside the time interval < at < tf . The last equality can be rewritten as 

d ^ W ^ = T J d d xdtth(x,at) (U(x,t) - P^(x,t) S je- l3W '^ , (A4) 

which vanishes because of equation (|4"Tj) . We thus have 

die-?™ 1 ) 
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Hence, the value of (e # Wj ) does not depend on a and can be evaluated by taking the limit a — >• which corresponds 
to an adiabatic evolution. But then it is well known, from classical thermodynamics, that Wj = —AF. This implies 
( e -0Wj) = exp(-AF). 



Appendix B: Supersymmetric formalism 

In this appendix, we use the superfield formalism, as explained in [17| and [29| , to rewrite the dynamical action and 
to interpret the invariances under the transformations (|50[) and (|52p in more compact and elegant language. In this 
formalism, the origin of these symmetries will appear more clearly. Besides, the effect of adding the Jarzynski term 
to make the time-dependent action invariant under (|52p . will also become more transparent. 

We introduce two anti-commuting coordinates 9 and 9 and define the superfield 

#(x, t, 9, 6) = 4>(x, i) + 9c(x, t) + c(x, t)6 + 66~4>{x, t) . (Bl) 
In terms of this superfield, the action <p 7 c, c), defined in (|49l) . can be written as 

£0, 4>, c, c) = J dMe £(#) with £(*) = T (5* D<f> + U{&)) , (B2) 

where the differential operators D and D are given by: 

D = (B3) 
86' K ' 

D = » _ £§° . (B4) 
89 r 8t ■ ' 

These two operators satisfy the anticommutation relations D 2 = D 2 = and {D, D} = — j^-§j ■ 
Integration with respect to the Grassmann variables is defined through the following rules [17| : 

J d9d9 1 = 0, J d9d9 9 = 0, J d9d9 9 = 0, J d9d9 96 = 1. (B5) 

(Integration and derivation are in fact identical). 

The action E is invariant under the transformations ([5U1) and (|52"|) which act by mixing the four fields (<j),4>,c,c). 
These two symmetries can be viewed as transformations of the superfield $ that leave the super-action £($) invariant. 

In the superspace formalism, the transformation (|50[) corresponds to an infinitesimal translation of the 9 coordinate, 
6 — > 6 + e. The generator of this transformation is given by 

<* = w (B6) 

Indeed, one can check that the superfield 5<f> = = 5(f>(x,t) + 95c{x 1 t) + Sc(x 7 t)9 + 965<j)(x,t) is given by 8<f> = 
ec(x, t) + e6<j)(x, t) = —c(x, t)e + 6<f>(x, t)e. If we identify each of the components we retrieve the transformation (|50p . 

Similarly, the transformation (|52[) corresponds to 9 — > 9 + 1 and t — ^ t + ^-e6. This transformation is generated by 
the operator 

The operators Q and Q that generate the supersymmetry transformations anticommute with D and D. Besides, 
they satisfy the anticommutation relations Q 2 = Q 2 = and {Q, Q} = y^-§i ■ When the potential U does not depend 
on time, the action £(*&) is symmetric under Q, and is invariant under Q upto a total derivative. This fact was 
checked in equations (|51I53[) and can be verified again using the supersymmetry formalism. If U depends explicitly 
on time the action £(<&) is not invariant anymore under Q. However, by adding to it the Jarzynski term (1261) . we 
obtain the modified action £j, defined in equation (|65p . which is invariant under Q. This property is manifest in the 
supersymmetric formalism in which the modified action is written as 

£j = T LD* + t + ^-99)) . (B8) 



r, 
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Appendix C: Supersymmetric Ward-Takahashi Identities 

When the invariances under the transformations generated by the operators (|B6|) and (|B7[) are implemented in 
the generating function Z(H,_H,L, L) defined in (j54")) . the Ward-Takahashi Identities (|55)l and (|5oT) are obtained. We 
follow closely the method of [24[ to derive these identities. In order to calculate correlation functions it is helpful to 
rewrite the sources as a superfield J, defined as: 

J = H + 9L + 9L + 99H (CI) 

where L(x,t) and L(x,t) are Grassmann fields. We thus have 

J d9d93{x, t, 0, 0)#(x, t, 9, 9) = H<j> + H4> + Lc + Lc. (C2) 

We note from this expression that H plays the role of an applied external 'magnetic' field. In this formalism, the 
generating function Z(H,H,L,L) becomes 

Z{3) = J V®eJ d*xdt«lfc»<-s:(*)+J*) _ (C3) 

(In the sequel, the integration element d xdtd8d6 will be omitted in general.) 

In order to derive the first Ward-Takahashi identity, we proceed as follows. In the functional integral (|C3I) . we 
make the change of variable $ — > $ with 3? = <I> + eQG>, where Q, defined in (|B6[) is the infinitesimal generator of 
the transformation (|50[) corresponding to 9 translations. Taking into account that the Jacobian is 1, we obtain 

Z(J) = J X>W- S W+ J * = J p$ e /-S(#+ £ Q#)+J(#+eQ*) = f p $e /-S(#+eQ#)+J(#+ £ Q#) ^ 

(The last equality results simply from the fact that $ is a dummy variable.) The fact that the action is invariant 
means precisely that J d9d9H(&) = J d9d9H(& + eQ<&). Therefore, we deduce that 

Z(J) = J V<f>eJ S+J * = J V^eJ -s+J*+J.Q* = J V ^ e I -E+J* (1 + e J jq$) . ( C5 ) 

This equation being true for any value of e we conclude that 

j V$> Ij & d xdtd9d9 JQ# J e$ ~ S + J * = . (C6) 

We now calculate explicitly the value of and substitute it in equation (|C6j) . 

/x>*{/d^^, t )a(«, t )-Z(,, t ) fet )}e/— ..o. (CT) 

By differentiating the generating function Z(3) with respect to H, we obtain 

J? {3 \ , = ! V*<j>(x a ,t a )expf-*+ J *. (C8) 
6H(x a ,t a ) J 

Similarly, we have 

52 - 52 a 62 X /no, 

-^c, ^c, and —^0. (09) 



Substituting these relations in equation (|C7I) allows us to derive the first Ward-Takahashi identity Jo 

^( B n- L w)- a - (C1 »» 

For the second invariance under the BRST2 transformation (|52p . we use the infinitesimal generator Q, defined 
in (|B7|) . After similar calculations, we find 



J V® jy d d xdtd9d93Q^ ef ~ S+J * = . 



(Oil) 
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After calculating explicitly JQG>, we obtain 



J £>* |y d d xdt Y~H( Xi t)c(x, t) + L(x, t)(<f>(x, t) - ^<j>(x, t)) + fT(x, t)c(x, t) | ^ - S+J * = . (C12) 

Expressing the fields in this equation as functional derivatives of the generating function Z, leads us to the second 
Ward-Takahashi identity ([SB]): 

d'sdi (Ih * 5 4 + L ( || - A * || U B S 4) = . (C13) 



r d*<5L V<5# T dt5Hj 5L 
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